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Abstract 

The Lifshitz formula is derived by making use of the spectral summation method which is a math- 
ematically rigorous simultaneous application of both the mode-by-mode summation technique and 
the scattering formalism. The contributions to the Casimir energy of electromagnetic excitations 



of different types (surface modes, waveguide modes, and photonic modes) are clearly retraced. A 



correct transition to imaginary frequencies is accomplished with allowance for all the peculiarities 
of the frequency equations and pertinent scattering data in the complex u plane, including, in 
particular, the cuts connecting the branch points and complex roots of the frequency equations 



(quasi-normal modes). The principal novelty of our approach is a special choice of appropriate 



passes in the contour integrals, which are used for transition to imaginary frequencies. As a result, 
the long standing problem of cuts in the complex oj plane is solved completely. Inconsistencies 
of some previous derivations of the Lifshitz formula are traced briefly. For completeness of the 
presentation, the necessary mathematical facts are also stated, namely, solution of the Maxwell 
equations for configurations under consideration, scattering formalism for parallel plane interfaces, 
determination of the frequency equation roots, and some others. 
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I. INTRODUCTION 



It is not overstatement to say that the Lifshitz formula is the basis for practically all the 
Casimir calculations dealing with the plane boundaries [lj. However, it was recognized in 
the literature {2-5] that the derivation of this formula in the original paper jf], Q| is very 
complicated. It proved that especially involved is the transition to the imaginary frequencies 
when obtaining the Lifshitz formula by deforming the contours of integration in the complex 
frequency plane. As far as we know, nobody has afterwards succeeded in repeating this 
part of Lifshitz' calculations. Furthermore, it was left unclear what are the excitations of 
electromagnetic field in matter which were taken into account in this formula. Now the latter 
point becomes important, in particular, in searching for the materials with given properties 
which are needed for obtaining the wanted characteristics of the Casimir forces. 

Later on, the Lifshitz formula was obtained by making use of different methods, namely: 



the quantum 



mation 



d theory technique in condensed matter physics 



16 



18], scattering formalism 



B3, 



mode-by-mode sum- 



15] , methods o 



19 



quantized surface modes and quantized Fresnel modes 

and so on. In the 



20], local Green functions [2l|, 



present paper we propose another derivation of the Lifshitz formula, namely, this formula 
will be obtained by making use of the spectral summation method. This approach is a 
mathematically rigorous simultaneous application of both the mode-by-mode summation 
technique and the scattering formalism. Inconsistencies of some previous derivations of the 
Lifshitz formula are traced briefly. 

Our derivation of the Lifshitz formula includes two steps: First we determine the spec- 
trum of electromagnetic excitations for given_boundaries taking into account the material 
properties of the media under consideration 



of the relevant zero point energies 



22 



24 



23] : afterwards, we accomplish the summation 



25] connected with all branches of this spectrum. 



For description of the electromagnetic fields in matter we shall use the dielectric formalism 
with allowance for the usual temporal dispersion. In the framework of this approach, one can 
concentrate attention on the electromagnetic field dynamics, while the dynamics of induced 
charges and currents is taking into account by introducing the permittivity and permeability 
e(u), /i(w). 1 



1 The electromagnetic field is coupled, through the Maxwell equations, with charges and currents; therefore, 
one can take, as the dynamic variable, the local displacement of a continuous charged liquid that describes 
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When determining the spectrum of electromagnetic field, all physically relevant solutions 
to Maxwell equations should be taken into account. Obviously, such solutions are all squared 
integrable solutions which belong to the L 2 functional space. In addition, this set of solu- 



tions should be complete. In rigorous spectral theory of differential operators [27|, [28|], the 
following assertion is proved: the functional space including the solutions which describe 
the bound states (discrete spectrum branch) and scattering states (continuous branch of the 
spectrum) is complete with respect to the L2-norm. It is these solutions that we shall find 
in the problem at hand and on this basis we shall determine the relevant spectrum, i. e. the 
admissible values of the electromagnetic oscillation frequencies u. 2 Upon determining the 
complete set of the solutions to the Maxwell equations, quantization of the electromagnetic 
field is trivial, because we are dealing here with the sum of two infinite sets (discrete and 
continuous) of noninteracting oscillators. The discrete set of oscillators describes the evanes- 
cent waves including the surface and waveguide modes 3 ; and a continuous set of oscillators 
is responsible for the scattering states. All these modes will be rigorously defined in the 
course of constructing the complete set of solutions to the Maxwell equations (see below). 

The summation of zero point energies of the discrete modes can be accomplished, ob- 
viously, in a straightforward way, while for the summation of such energies appertaining 



to scattering states one has to take advantage of the scattering formalism [3JJ. Preced- 
ing derivations of the Lifshitz formula in the framework of the mode-by-mode summation 
alone or by making use of the scattering formalism only are apparently inconsistent from 
the mathematical point of view. Strictly speaking, each of these methods can be applied 
separately to one branch of the spectrum, namely, the mode-by-mode summation is appli- 
cable to discrete part of the spectrum solely and the scattering formalism enables one to 
take into account the scattering states, i. e. it is applicable to the continuous part of the 
electromagnetic spectrum only. 

The understanding of this fact was absent in previous papers, for example, in Ref. [kJ 
applying the mode-by-mode summation technique it was noted that "... the surface modes 



the free electrons inside the medium [26 

2 In the general case, the spectrum of a differential operator may be more complicated [29j . However in the 
Casimir calculations such problems do not occur. 

3 In physical problems dealing with the plane boundaries the surface modes and waveguide modes are called 
sometimes the guided modes [3o| . 
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alone do not yield right expression" for the correct retarded dispersion force. But the author 
didn't find out what should be added to the surface modes in order to obtain correct result. 
He claimed further that the calculations presented in papers 



11 



12| are incorrect because 



the cuts in the complex frequency plane u were not taken into account when applying the 
argument principle [3, [ijj]. For solving the problem of branch cuts, it was proposed to 



consider the plates of a finite thickness instead of semi-infinite slabs {3 , 4], |l3j, [14 ] . 

In fact in papers jll[ Q and in many others [l[ I13I-I8]. besides neglecting the cuts in 
the complex frequency plane, that is a mathematical mistake, there is another physical 
inconsistency, namely the continuous spectrum was completely ignored. The author of Ref. 
12[ argued erroneously that ". . . the well-known Lifshitz formula for the retarded van der 



Waals attraction between two solid half-spaces can be obtained from the zero-point energy 
of the interacting surface plasmons." Keeping in mind the well known paper 10], one 



could make wrong inference that both nonretarded and retarded van der Waals interaction 
given by Lifshitz formula are due to nonretarded and retarded surface plasmon interactions 

n 

respectively [32| . Thus it turned out as if the electromagnetic bulk excitations (scattering 
states) are irrelevant to the calculation of the vacuum forces. 

On the other hand, when deriving the Lifshitz formula in the framework of the scattering 
formalism alone (see, for example, [331]). an artificial technique was applied in order to take 
into account the evanescent waves (bound states), namely, analytical continuation of the 
reflection and transition amplitudes to the corresponding region |34j. 

As a matter of fact, it will be shown in the present paper that all these problems can be 
solved consistently only by a simultaneous allowance for the contributions to the vacuum 
energy generated by the evanescent modes (discrete branch of the electromagnetic spectrum) 
and by the scattering states (continuous part of this spectrum), provided the appropriate 
mathematical methods will be employed for treating the different parts of the spectrum. 

The important step in the Casimir calculations is the transition to the imaginary frequen- 
cies. In our opinion this transition was not rigourously justified at least in the framework 
of the mode-by-mode summation technique. It is this issue that resulted in foregoing con- 



fusion and erroneous c 



summation ll|, |3j, |4 , 



QflQ 



10 



aims in the first attempts to derive the Lifshitz formula by the mode 



181 ] . In order to accomplish this transition in a consistent way the 



analysis of the analytical properties of the frequency equations and the scattering data in 
the complex u plane should be conducted preliminary and the needed cuts connecting the 
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branch points in this plane should be done. In the pertinent literature [l-4, llOM22j| these 
issues were not examined. We are going to fill in this gap. 

It is worth noting here that the analytic properties of the scattering matrix (or the Jost 
function) in the Casimir studies are different in comparison with the standard theory of 
potential scattering. In fact they are close to those for the Klein-Gordon equation, the 
role of mass squared being played by k 2 , where k is the wave vector along the unbounded 
dimensions. This implies, in particular, that the analytic properties of the scattering matrix 
with respect to the complex frequency u in the Casimir calculations should be revealed by 
direct analysis of its explicit form without referring to the nonrelativistic potential scattering 
in the framework of the Schrodinger equation. 

The layout of the paper is the following. Section |TT] is devoted to the spectrum of elec- 
tromagnetic excitations in the problem under consideration. For the Casimir studies the 
materials permitting the surface waves are of current interest. The permittivity of these 
media can acquire, in certain frequency bands, negative values. Typical examples of these 
materials (metals and isotropic dielectrics) are discussed briefly. In Sec. Ill Al the general 
solutions to the Maxwell equations for plane parallel interfaces are constructed. The both 
branches of the spectrum are considered in detail. In Sec. II III the Lifshitz formula for the 
Casimir energy is derived by summing up the zero point energies of discrete modes and 
scattering states. In Sec. IIVI (Conclusion) the obtained results are summed and unsolved 
problems in this field are briefly outlined. 



II. THE SPECTRUM OF ELECTROMAGNETIC EXCITATIONS 

The electromagnetic excitations in the media are diverse, and they essentially depend on 
the dielectric and magnetic properties of the background. 

It is worth mentioning that we are interested in the long wave excitations with A ^ d, 
where d is the characteristic scale in the Casimir calculation (c? ~ 10 - 100 nm) and A is the 
considered wave length. The description of the condensed matter in terms of the dielectric 
permittivity assumes the long wave approximation as well. 

In applications the dielectric properties of the materials play the leading role, therefore 
in what follows we put the magnetic permeability equal to a unit. 

Of particular interest for theoretical and experimental Casimir studies are the media with 
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dielectric permittivity taking on negative values over a certain frequency range. Only inter- 
faces of such media support surface electromagnetic waves which can be treated as collective 
excitations of electron density. The surface waves together with waveguide solutions belong 
to the discrete branch of the spectrum (see Sec. Ill Bl) . 

Metals are typical media of this sort. The overall picture of the electromagnetic waves 
(excitations) in metal resembles that in plasma but is not identical with it. In the bulk of 
the metal the local electron density oscillates with the characteristic plasma frequency 

u 2 p = , 2.1 

where n is the mean electron density in the metal and m is the electron mass, the role of the 
restoring force being played by the electrostatic field. The volume plasma excitations give 
rise to surface excitations of the electron density (surface plasmons) which propagate along 
the metal-dielectric or metal- vacuum interface. Obviously, the oscillations of the free charge 
density cause, according to the Maxwell equations, the oscillations of the electromagnetic 
field. 

For metals, the following dielectric function is commonly used (plasma model) 

e(w) = l-^, (2.2) 

where u p is the plasma frequency (12.11) . The dielectric permittivity e(u) takes negative values 
in the region u < oo p . The field oscillations caused by the surface oscillations of the electron 
density are localized near the surface as well. They are sometimes called surface plasmons 
(see the corresponding analysis of the solutions to the Maxwell equations in Section Til A\i . 

Inside the metal the bulk waves propagate with the frequencies u > u p . The transversal 
and longitudinal bulk waves form the continuous branch of the spectrum. 

In an isotropic dielectric the permittivity is described, with a good accuracy, by an oscil- 
lator model |35] . To gain better understanding, one can consider the atoms in a dielectric as 
weakly damped harmonic oscillators with eigenfrequencies Ui excited by an external electric 
field E = Eoe~ iW *. The iV-oscillator model gives the following dielectric function: 

N S u 2 

«H =«(«>) + £ J, i (2.3) 
j J J 

where Sj is the j-th oscillator strength and Fj is its relaxation parameter. In what fol- 
lows we neglect the absorption (undamped oscillators, Tj = 0) and consider real dielectric 
permittivity. It is negative for u approaching ujj from the right. 



Thus, the condition for the existence of the surface plasmon at the flat dielectric-vacuum 
interface is fulfilled near the narrow absorption lines of the dielectric media. The number of 
the plasmons in this case is equal to the number of oscillators in the model. 

The bulk waves in dielectrics (continuous spectrum) lie in the ranges of frequencies where 
e((j) is positive. 



A. General solutions to Maxwell equations 

First we recall briefly the formulation of the Maxwell theory for compound media with 
constant permittivity e and permeability \i in each region. In this case, the harmonic in 



time electric (E) and magnetic (H) fields are described by the Maxwell equations [36] 



VxE = i-B, V ■ D = 0, (2.4) 

c 

VxH = -i-D, V-B = 0, (2.5) 

c 

D = eE, B = /iH, x^S, (2.6) 

which hold outside the interface £ separating the regions with different e and /i, and by 
matching conditions on S. These conditions require the continuity of tangential components 
of the fields E and H when crossing £ 

discont (E||) = 0, discont (Hy) = . (2.7) 

The Gauss units are used and it is assumed that external charges and currents are absent 
(both volume and surface ones). The common time factor e~ luJt will be dropped. 

When allowing for the time dispersion, the permittivity e and permeability /i in material 
equations (12. 6j) and further should be treated as the functions of the frequency u. 

General solution to Eqs. (12.41) . i \2.5L and ( 12 .61) can be represented in terms of two inde- 



pendent Hertz vectors in the following way 



36 



E = V x V x IT, H = -is - V x IT (TM-modes); (2.8) 

c 

E = in - V x n", H = V x V x II" (TE-modes) . (2.9) 

c 

Here II' is the electric Hertz vector, II" is the magnetic Hertz vector, and c is the velocity 
of light in vacuum. 
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In each region with given e{ui) and /u(w) the Hertz vectors obey the Helmholtz vector 
equation 

(v 2 + P)n = o, (2.10) 

where the wave number k is given by 

tu 2 

k 2 = e(u)fji(u) -r . (2.11) 

The boundary conditions (12. 7p are given for the complete fields E and H, i.e. for the 
sums of the TE- and TM modes. However, for the interfaces of high symmetry the TE- and 
TM-modes, fortunately, decouple in (12. 7p . 

It is known that in the source-free case the general solution of Maxwell's equations can be 



derived from two real scalar functions 



39 



- |4lj which may be chosen in deferent ways. Thus 



the Hertz potentials have in fact only one nonvanishing component. The precise choice of 
these components is determined by the interface geometry and the reasons of simplicity and 
convenience. 

We consider the flat interface parallel to the (x, y) coordinate plane and cutting the z-axis 
at the point z = zq. Let e^, e y , e z be the unit base vectors in the chosen coordinate system. 
The ^-components of the Hertz potentials are treated as two independent functions in the 
general solution to the Maxwell equations 

IT = e,e iks $(z), II" = e,e iks ^(^) . (2.12) 

Here k is a two-component wave vector parallel to the interface, k = (k x , k y ), and s = (x, y). 
The common time-dependent factor e~ tljjt is dropped as usual. 

Substituting the Hertz potentials ( 12TT21 in §ZM) and O we get the fields for the TE- 
modes 

E = fi^-k^y^^iz), (2.13) 
H = ifc x e a .e iks ^'(z) +e z k 2 e iks ^{z) (2.14) 

and for the TM-modes 

E = ik x e x ^ s &(z) + e 2 fcV ks $(z) (2.15) 
H = -e^k x e y (^ s §{z) . (2.16) 



S 



Without loss of generality, we directed the x-axis along the vector k: k = (k x = ±k, 0). 

For the fields ( 12TT3]) -( 12~T6]) to obey the Maxwell equations the Hertz potentials 

should meet the Helmholtz equation (I2.1(jp which now assumes the form 

- <$>"{z) = (efi^ - k 2 ^j $(z) (TM-modes), (2.17) 

-*"(*) = (en^ -k*)^(z) (TE-modes), (2.18) 
—00 < z < 00, z ^ z , k 2 = k 2 + ky . 

Substitution of the fields (I2.13I) - (I2.16I) into the continuity equations (12 .7p results in the 
matching conditions at z = z separately for the functions $ and ^ 

[e(z )$(z )) = 0, [&(z )} = (TM-modes) , (2.19) 
Mzo)tt(zo)] = 0, [*'(zo)]=0 (TE-modes), (2.20) 

where the notation 

[F(z)]=F(z + 0) - F(z-O) (2.21) 

is introduced. 

For a given value of k 2 > the differential equations (12.17j) . (12. 18ft with matching con- 
ditions (I2.19p . (I2.20p and physical conditions at infinity z —> ±00 result in two spectral 
problems for the TE and TM polarizations. The frequency 00 plays part of the spectral 
parameter. For complicated functions e(u) and the spectral parameter may enter into 
the initial differential equations nonlinearly. 



B. The branches of electromagnetic spectrum relevant for calculation of vacuum 
energy 

In order to calculate the vacuum energy of the electromagnetic field, one has preliminarily 
to quantize this field with allowance for given boundary conditions. To this end, the full set 
of solutions to the Maxwell equations including all physically relevant ones is needed. In the 
directions parallel to the interfaces the electromagnetic field dynamics is free. Therefore, 
we are concerned with the behaviour of the solutions normal to the boundaries, along the z 
axis. 
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As was mentioned in the Introduction, the full set of physically relevant solutions in the 
problem under consideration comprises discrete natural modes and scattering states. 

Normal modes are the solutions to the Maxwell equations which are localized near and 
between the interfaces, their energy being localized too. They are analogous to the bound 
states in quantum mechanics. These solutions are square-integrable on the whole z axis. 
The corresponding frequencies (the eigenvalues of the spectral problem) take on discrete real 
values. 

The scattering states (propagating modes) are described by incident waves and outgoing 
scattered (or reflected and transmitted) waves. Their frequencies are real and positive 
forming a continuous spectrum. 



1. Discrete part of the spectrum 

We start with considering a single flat interface in order to elucidate the properties of 
dielectric and magnetic functions required for the existence of normal modes in the problem 
under study. 

a. Normal modes for a single plane interface. Let the plane z = separate two uniform 
half-spaces characterized by ei, \i\ and e 2 , or symbolically (ei, (J,i\e2, ^2)- The general 
solution to the Maxwell equations are defined by two functions <&(z) and ^f(z) obeying (12. 17ft . 
( 12.18!) an d the matching conditions ( 12.191) . ( I2.20p . The eigenfunctions in the problem should 
decrease exponentially in both directions from the interface, i.e. they behave like outgoing 
waves with imaginary wave-vectors 



A ie ~ iklZ , z < 0, 
= < (2-22) 

A 2 e ik * z , z>0, 

where A\ and A 2 are the constant amplitudes, 

22 2 
k 2 a = e a fx a ^ - k 2 = ^ - k 2 < 0, c 2 a = —, a = 1,2, (2.23) 

C C a E a [l a 

provided that 

k a = +i I k a I , a = 1,2. (2.24) 

On substituting (12.221) and the analogous representation for the function ^(z) into the 
matching conditions fl 2 . 1 9 [) and ( I2.20p we arrive at the following equations determining the 
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eigenfrequencies 



e 1 k 2 + e 2 k 1 = (TM-modes), (2.25) 
H^+^h = (TE-modes). (2.26) 

In view of condition (I2.24p the frequency equations (12.251) and (12.261) may have real roots 
(real eigenfrequencies) only if the permittivities £i, e 2 and permeabilities fJ>i,^ have different 
signs. 

This condition is fulfilled, for instance, at the metal-vacuum interface, provided the di- 
electric function of the metal Ei(u) is described by the plasma model (12.21) . In this case we 
put Ei — ii\ — Hi — \. Substituting all this into (I2.23j) and (I2.25p . one obtains the dispersion 
equation for the surface plasmon in the TM-polarization 

fi 4 -(l + 2g 2 )fi 2 + g 2 = 0, (2.27) 

where the dimensionless variables 

0=-, g=£, k p = ^ (2.28) 

are introduced. Solution to Eq. (I2.27P is given by 




V(q) = \ h + l'-Xh + l'- ( 2 - 29 ) 



The dispersion curve ft(q) is presented in Fig. [Ik, Another solution to Eq. ( I2.27P leads 
to positive k\ that corresponds to propagating waves in medium 1 (dotted curve in the 
upper-left corner in Fig. [T^l). In view of the asymptotic behavior 

n{q)^q(l-j + --\ q^O, Q{q) -+ ±= (l _ J_ + . . , q oo , ( 2 .30) 

the dispersion curve Vt(q) tends in these limits to two straight lines Vt = q and Q = 1/V2 
from below. Here VL = 1/ v^2 is the dimensionless frequency of the surface plasma oscillations 
(see Fig. [Tk). 

With /ii = /i2 = 1, obviously, there is no surface modes in the TE-polarization. 
The surface waves can be sustained by the dielectric-vacuum interface too, as the dielec- 
tric permittivity, with account for the dispersion (12. 3p . may take on negative values. The 



corresponding dispersion law can be found, for example, in 

11 



42). 



b. Normal modes for two parallel plane interfaces Now let us consider two interfaces, 
parallel to the xy plane and cutting the axis z at the points z = —a and z = a. For simplicity 
we consider at first symmetric configuration that can be symbolically represented as 

, \l\ | E 2 , /i2 I £1, A*i) • 

For given polarization, the normal modes can be subdivided into symmetric or asymmetric 
solutions. Thus, we have for the TM-modes 



Ae- ikl( - z+a \ z<a; 
±Ae ik ^ z ~ a \ z>a, 



and the analogous ansatz for the function *ff(z) in the case of the TE- modes. The wave 
vectors k\ and k 2 in (I2.31[) are defined, as before, by Eq. ( 12.231) . For \z\ > a we again 
consider the 'outgoing waves' with the imaginary wave-vector k\ < 0, while in the region 
\z\ < a the wave vector k 2 may be both imaginary, k\ < 0, or real, k\ > 0. In the first 
case, the waves are just surface ones localized near the interfaces and in the second case, 
they are the waveguide solutions 4 describing the standing waves between the interfaces and 
the evanescent waves outside this region. Indeed, both solutions are normal modes because 
they are squared integrable on the whole axis z. 

The matching conditions at the points z = —a and z = a result in the respective frequency 
equations. For the TM-modes these equations read 

L a = ikis 2 + £1^2 tan(aA;2) = (asymmetric modes), (2.32) 
L s = ik\E2 — £ik 2 cot{ak 2 ) = (symmetric modes). (2.33) 

When defining the spatial symmetry of the electromagnetic modes we follow the papers 
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45], namely, the symmetric TM-modes are constructed by making use of the symmetric 
function <&'(z) and, consequently, antisymmetric function $(2) in Eqs. f)2.15p . f)2.16p . It 
implies that the longitudinal component of the electric field E x (z) is symmetric with re- 
spect to the plane z = and the transverse components of the fields E z (z) and H y (z) are 
antisymmetric. The antisymmetric TM-modes have the contrary symmetry properties, i. e. 



In 



331 ] . such solutions are referred to as the cavity modes. 
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E x (z) is antisymmetric with respect to the z = plane, and E z (z) and H y (z) are symmetric. 
In the literature HQ reverse notation is used too. The eigenfrequencies of the symmet- 
ric modes are denoted by u s (k) and those for the antisymmetric modes by u a (k). In the 
literature different notation for these frequencies is used, for example, sometimes u + (k) is 
referred to the higher frequency modes and u~(k) applies to the lower frequency collective 



oscillations 



43 



M\- 



By making use of the trigonometric relations 

2 tanx 



tan2x 



tan 2 x 



2 cot x 
cot 2 x 



1 cot x — 1 

the frequency equations (12.321) and (I2.33j) can be combined into one equation 

2^7 e 2 ki 



tan(2afc 2 ) 
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1 + 7 2 ' ' Eik 2 

The frequency equations (12.321) and (I2.33P can also be rewritten in the exponential form 



(2.34) 



1 + r u e 



2iak2 



(L a = 0), 



rue 



2iak2 



(L, 



(2.35) 

0). (2.36) 

Multiplying together these two equations we arrive at the exponential form for the total 
frequency equation (12.3411 

D(u) = 1 - r 2 2 e 4Mfc2 = (L a L s = 0), (2.37) 



where the notation 



r 12 = £A ~ £l f 2 = 2^1 (2.38) 
e 2 k\ + E\k 2 7 + 1 

is introduced. It will be shown later (see Eq. (I2.50p ). that r 12 is the reflection amplitude for 
a single interface. 

The frequency equations for the TE-polarization is obtained from the above equations by 
replacement s a — > fi a , a = 1,2. 

To comprehend the structure of the discrete spectrum in the presence of two interfaces 
we address typical examples. The solutions to the frequency equations (12.321) and (12.331) are 



studied in detail in the theory of dielectric waveguides 
in the theory of surface plasmons 



42 



44 



50 



46[ and optical waveguides [48l. |49| . 



52| and so on. However, in these applications 



the roots of the frequency equations are often sought of in the form k = k(u>) (wave vector 
as function of frequency). We are interested in the inverse function oj = u(k). 
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While studying the frequency equations (I2.32p and (I2.33|) their plots in terms of the 
variables k\ and fc 2 instead of the initial variables k and u are helpful and transparent. To 
this end, the considered frequency equation is appended with the relation between k\ and 
k\ following from their definition (I2.23p . At fixed u one obtains two equations connecting 
k\ and & 2 - The intersection points of the curves defined by the left-hand sides of these 
equations give the solutions of the initial frequency equation (see, for example, [46l. |53|). 

Let us briefly describe the structure of the discrete spectrum of electromagnetic exci- 
tations for some configurations. Not all these examples can have a direct relation to the 
experimental studies of the Casimir forces. Nevertheless, all of them turn out to be in- 
structive in revealing the general properties of the electromagnetic spectrum for two plane 
parallel interfaces. 

We specify the configuration I in the following way: 

£1 = A*i = 1, £2 and /i 2 are some constants (2.39) 
(a material plate in vacuum without dispersion). This configuration is studied in detail in 



radio engineering as the simplest example of a waveguide 37], |46j]. When the conditions 



s 2 > 0, n 2 > 0, £ 2 /i2 > 1 (2.40) 

or 

e 2 < 0, n 2 < 0, e 2 /i 2 > 1 (2.41) 

are satisfied, there exist waveguide solutions in the form (12.31!) with k\ < and k\ > for 
both the TE- and TM-polarizations. For a given value of k 2 , the number of such solutions 
is finite and is increasing with u. The respective real eigenfrequencies lie in the interval 
(c 2 k, Cik) with ci = c. 

As expected, the surface modes (k\ < 0, k\ < 0) in the present configuration exist only 
under the conditions (I2.4ip . in other words, if the plate is made of the left-handed material 



(both 62 and // 2 are negative 54J, |55[). In the general case, the number of modes of this kind 
is equal to 4: there are 2 surface modes with the TE-polarization (symmetric oscillation and 
antisymmetric oscillation) and 2 analogous modes with the TM-polarization. 

For symmetric surface modes it is enough to fulfil the condition (I2.4ip . For antisymmetric 
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surface modes, one has to impose, in addition to (I2.4ip . the following conditions 



( ,..- I \e2\ 1 TM-polarization; 

< — (e 2 // 2 - 1) < I (2.42) 

' | /X2 1 1 TE-polarization. 

The real frequencies of the surface modes lie in the interval (0, 02k). 

As the next example we consider configuration II which is complimentary to the pre- 
vious one, namely, £1 and /ii are some positive or negative constants such that £i/ii > 1 
and £2 = ^2 = 1 (two semi-spaces of the same material without dispersion separated by the 
vacuum gap). 

This configuration is not investigated in radio engineering and wave guide optics but it 
is closely related to the Casimir calculations. In the case of left-handed material, obeying 
the conditions 

£1 < 0, /ii < 0, £i/ii > 1 , (2.43) 

there exist surface symmetric and antisymmetric waves in both polarizations 
one of the branches of the hyperbola defined by the equation 

,2, ,2 



55, 



561 ] . Indeed, 



a k 2 — a k, = — — (£i//i — 1) > , k a — in a , a = 1, 2 (2.44) 

cr 



at negative £1 and \i\ always crosses in the first quadrant the curves defined by the equations 

CbKi = — — (an 2 ) tanh(a/«2) , cik>i = — — (0^2) coth(a/« 2 ) 
^2 £2 

and two more obtained by replacing e a by fi a , a = 1,2. Therefore, the discrete spectrum 

for the present configuration exists only in the case of left-handed materials and consists of 

4 surface modes. The frequencies of these modes belong to the interval < u < c±k. 

Keeping in mind the possibility of total reflection at the interface of the media 1 and 2 
one can try to find the solutions for this configuration with k\ > and k\ < 0. However 
it is easy to see that there are no such solutions with real frequencies as the equations are 
complex. Even if such solutions existed, they are not eigenfunctions of the discrete spectrum 
because they are not squared integrable on the infinite axis z. 

For the configuration under consideration there are no waveguide solutions to the Maxwell 
equations. Indeed, for the existence of such solutions it is required that k\ < 0, and k\ > 0. 
From these conditions we arrive at the conflicting inequalities < u> < c\k and u> > ck, with 

Cl = C/ y/Eifll < C. 
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Now let us consider the examples of the discrete electromagnetic spectrum explicitly 
allowing for the temporal dispersion in the media. As configuration III we take a material 
(metal) slab in vacuum, the dielectric properties of the slab being described by the simple 
plasma model (12. 2p (see Fig. [Tb): 

e 1 = fx 1 = l, £2 {u) = l-% /i 2 = l. (2.45) 
This configuration supports two surface waves in the TM polarization. These modes are 



well studied theoretically 
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55| and experimentally, for instance, in measurements of the 



energy losses by the electrons passing through the metal films 43|, |44J. There is a vast 



literature on this subject 



42 



50 



52|. Our citing is minimal. 
The physical picture of these collective excitations is the following. At both boundaries 
of the metal plate (z = ±a) with a — > oo there exist plasma oscillations of the electron 
density with the same frequency u p /y2 and dispersion law (12.291) . If the thickness of the 
plate or film is finite these oscillations interact reciprocally and the energy level (I2.29P splits, 
as usual, in two levels (asymmetric and symmetric ones) with frequencies uj a {k) and u s (k) 
defined by Eqs. (I2.32|) and f!2.33j) . respectively (see Fig. [Tb). 

It is important that both curves u s (k) and u a (k) are placed to the right of the line u = ck 
and below the horizontal line u = u p . Thus the eigenfrequencies of the surface modes lie in 
the interval < oo s (k), u a (k) < ck for < k < oo. 

The curve u a (k) at all k lies above the curve u s (k). Near the origin both the curves 
touch the line u = ck from the right. With k — > oo the curves u s (k) and u a (k) tend to 
the horizontal line u = u p / a/2, respectively, from below and from above. For 2a — > oo the 
frequencies u s (k) and u a (k) tend to the universal curve (I2.29P from below and from above. 

The characteristic value of the longitudinal wave vector k in the problem under consider- 
ation is k p = uj p /c. It is within the vicinity of k ~ k p where the "bend" of the curves u s (k) 
and oj a (k) is observed. 

It is evident that there are no waveguide solutions (kf < 0, k\ > 0) in the present case. 
Indeed, from the definition of the wave vectors k\ and k 2 (see Eqs. ( 12.231) and (12.21) ) we 
obtain 

— k\ + k\ = ^ . 

cr 

For the waveguide solutions the left-hand side of this equation should be positive. Hence, 
this equality cannot be satisfied. 
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FIG. 1: (color online) Dispersion curves for the discrete modes in the TM-polarization for different 
configurations: single metal-vacuum interface (a); metal plate of thickness L = 2ak p = 0.63 in 
vacuum (b); vacuum gap of the same width L = 0.63 in the metal bulk (c); vacuum gap of the 
width L = 6.3 in the metal bulk (ci). Frequencies Q, to, and the wave number q are presented in 
the dimensionless units defined in Eq. (I2.28D . The gray filled regions correspond to the continuous 
spectrum. 



Now we consider configuration IV, complimentary to the previous one, namely, a vac- 
uum slot, 2a in width, in the bulk metal described by the plasma model (see Fig. [Tfc,d) 

,2 



£i = 1 |, /^l = 1, £2 = A*2 



1. 



(2.46) 



First, we consider the surface waves in the present set-up, k\ < 0, k\ < 0, which exist 



in the TM polarization 



451 ] . Again, we have two such collective modes, symmetric and 



antisymmetric ones, their frequencies being placed in the interval (0, ck). The lay-out of the 
dispersion curves is interchanged in comparison with the previous case, oo s (q) > u a (q). The 
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curve u) a (q) lies below the horizontal line u = l/v2 and to the right from the line u = q. 



At the point q = q CI = l/yl + L/2 the dispersion curve u s (k) crosses the light line u = q 



and becomes the lowest waveguide mode (see Fig. [Tb,d) and 45j|). 5 Thus the symmetric 
surface mode oo s (q) becomes in this configuration a 'hybrid' one. It is this fact that leads to 
the known result, namely, the contribution of the surface modes to the Casimir energy for a 
given configuration proves to be considerable only at small distances. Indeed, because of the 
hybrid nature of the symmetric mode, the dispersion curves u s (q) and uj a (q) move apart in 
the region q < q cr (Fig. [Tb,d), this divergence being maximum at small L (compare Fig. [lb 
and Fig. [TH). Obviously this implies in turn that the mutual interaction of these collective 
excitations is maximum at small L. 

With increasing the gap width L, in this configuration there appear waveguide solutions 
{k\ < 0, k\ > 0) satisfying the conditions 

c 2 k 2 <u 2 wg (k) <u 2 p + c 2 k 2 . (2.47) 



Thus, they lie to the left from the line u = ck and below the curve u = \foj p + c 2 k 2 . 



Approximate 



are given by 32] 



y the frequency of these modes and their total numbers, for each polarization, 



u 2 (k) ns k 2 c 2 + n 2 7T 2 c 2 /4a 2 } n = 0, 1, 2, . . . , n < 2au p /nc . (2.48) 

These modes exist only at large gaps (compare Fig. [lb and [Thl) . Indeed, for gold u p = 
1.38 ■ 10 16 s" 1 . Therefore, the dimensionless width of the gap L = 0.63 implies 2a = 14 nm 
and L = 6.3 corresponds to 2a = 140 nm. In the first case, there are no pure waveguide 
modes (Fig. [lb) and, at the same time, the surface modes interact strongly because their 
dispersion curves are distinguished considerably. In the second case (Fig. [TH), only two pure 
waveguide modes appear while the surface modes at q > 1 practically coincide with each 
other and with the single interface dispersion curve fl(q). The latter means that in this region 
the surface modes do not couple. For copper, we have close figures: u p = 1.46 ■ 10 16 s -1 , 
2a = 12.6 nm (L = 0.63), and 2a = 126 nm (L = 6.3). 

The fact that waveguide modes exist only for large gaps suggests that they contribute 



appreciably to the dispersion force at large distances 32|, where contributions of the surface 



Probably, this implies that the symmetric surface TM-mode in this configuration may directly be excited 
by the light beam contrary to the plate configuration, Fig. [Tb [42 - 44 , 5Cj . 
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modes practically vanish. 

The examples studied above reveal general properties of the discrete spectrum in the 
problem at hand and permit us to establish the lower boundary of the continuous spectrum. 
Let us summarize these properties: 

i) The discrete spectrum may comprise the solutions of two types: the surface modes 
and the wave guide modes; 

ii) At fixed k the frequencies of the surface modes (if they exist) lie in the interval 
(0, u-(k)), while the frequencies of the waveguide solutions (if they exist) belong to the 
interval (u-(k),u + (k)), where the boundary values uj-(k) and ou + (k) (u-(k) < ou + (k)) 
are defined by the explicit form of the dielectric function e(u); 

iii) The lower boundary of the continuous spectrum is defined by the condition k\{ui) > 0. 
Here we list the values of u + (k) and oj-(k) for the configurations I - IV. 

civ 

I: ou-(k) = c 2 k = , u)+(k) = ck, c 2 < c . 

In this case the discrete spectrum may include the surface modes and the waveguide modes 
with the frequencies from the above specified intervals; the continuous spectrum u is defined 
by uj > oj + {k). 

II: uj-(k) = c\k = , u)+(k) = ck, c\ < c . 

The discrete spectrum may consist of the surface modes only with the frequencies from the 
above specified interval; the continuous spectrum u exists for u > u-(k) = c±k. 

HI: w_(Jfe) = ck, cu + (k) = ^c 2 k 2 + uj 2 p . 

The discrete spectrum is formed by the surface modes with frequencies from the interval 
specified above; the continuous spectrum is defined by oj > u-(k) = ck. 

IV: u^{k)=ck, u+{k) = ^Jc 2 k 2 + u 2 . 

The discrete spectrum comprises the surface modes and the waveguide modes with fre- 
quencies in the corresponding intervals; the continuous spectrum exists for u > ou + (k) = 
sjc 2 k 2 + uj 2 . 
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In what follows it is significant that the points U-{k) and uj + (k) are the branch points 
for the root functions k\{u) and k 2 {u) (see Sec. IHI Al) . 

The configurations considered above are unbounded, i.e. open. As a result in addition 
to the real roots the frequency equations have also discrete complex roots u(k) = u'(k) + 
iu"(k) which correspond to the quasi-normal modes 3l|. These modes are not squared 
integrable and, therefore, they are not normalizable. That it why the quasi-normal modes 
are not included into the completeness relation. However their existence should be taken 
into account when choosing the contours for integration in the complex u plane (see Sec. 

MM- 

It is worth noting here that the quasi-normal modes are explored in fact in the experiments 
dealing with the radiating plasmons 44|. Usually, in this field the complex wavenumbers 
are considered k(u) = k'(u) + ik"(u) instead of the complex eigenfrequencies. 



2. Treatment of the continuous part of the spectrum 

To study the contribution of the continuous spectrum to the vacuum energy, the scattering 
matrix is required, i.e. the S-matrix for the wave equations (I2.17P and (12.181) with the 
matching conditions (I2.19p and (I2.20p . respectively. Here we deal with one-dimensional 



scattering on the infinite z-axis, — oo < z < oo. This problem has some peculiarities [57H60]. 
Specifically, one has to consider the incidence of the initial waves from the left and from the 
right separately. Thus, we must add such initial waves to the outgoing and standing waves 
in Eqs. (12.221) and (I2.3ip considered in the pertinent eigenvalue problems. The respective 
reflection amplitudes, r;, r r , and transition amplitudes, ti, t r , make up a (2 x 2) matrix that 
plays the role of the scattering matrix in the problem under consideration 



S{u) = (2.49) 




a. Single plane interface In the case of the one plane interface (ei, fii | s 2 , [i 2 ) we have 
for the TM-polarization 

^2^1 — ^1^2 / \ 

n = r n = — — ■ — , r r = r 21 = -r X2 = -r t , (2.50) 

tl = tl2 = e 1 h+ 1 e 2 k 1 i U = hl = ej'tl^ ^^^i- (2.51) 
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The analogous amplitudes for the TE-polarization are obtained from Eqs. (I2.50P and (I2.5ip 
by the substitution e a — » fi a , a = 1,2. 

The amplitudes t and r for the TM- and TE- polarizations have the poles at the points 
where the variable u coincides with the eigenfrequencies defined by Eqs. (I2.25P and (I2.26p . 
This property is preserved for several interfaces too (see below). The single plane interface 
is non symmetric configuration, therefore the ^-matrix (I2.49P in this case does not possess 
appealing physical properties. 

b. Electromagnetic scattering by two parallel plane interfaces Now we address the scat- 
tering matrix for two plane interfaces which are parallel to the xy plane and cut the z 
axis at the points z = —a and z = a, in other words, we envisage the configuration 



(ei,/ii | e 2 ,fJ, 2 I £3,^3)- The exp 
conditions fl2TT9|) and f l2~20]) (35, 



icit form of t and r is again derived from the matching 
Taking careful account of the phase factors at the 



points z = —a and z = a one obtains 

ri2 + r 23 e iik2a _ t 12 t 23 e 2ik2a 
n = ri23 = 1 + r 12 r B ^ ' tl = h23 = 1 + r 12 r 23 e^ ' ^ 

r 32 + r 21 e 4ik2a _ t 32 t 2l e 2 ^ a 

Tr ~ 1 + r 32 r 21 e 4 ^ ' r ~ 321 1 + r 32 r 21 e 4 ^« ' 1 ' j 

where r^ and t^ are defined in Eqs. f)2.50p and (12.511) . 

Further we consider the special case of (12.521) and f ]2.53[) . the symmetric one, (ei,fi\ \ 
£2,^2 I In the symmetric set-up one has 

_ r 12 -r 12 e ti 2 t 2 ie , . 

r = ri = r r = r\ 2 \ = 7, — n — , t = ti = t r = tioi = 7, — 7-, — . (2.54) 

The amplitudes r and £ in Eq. (I2.54p have the poles at the points of the discrete spectrum 
(see Eq. (12.371) ). and they obey the relations 

\t\ 2 + \r\ 2 = 1, (2.55) 
ri+ft = 0, (2.56) 

where the bar means complex conjugation. Now the scattering matrix ( I2.49P is given by 

S(u) = K T \ . (2.57) 

Due to Eqs. (12. 55ft and (12. 56ft the 5-matrix ( 12. 57ft is a symmetric unitary matrix 

SS* = S*S=1. (2.58) 
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In the case of potential scattering on an infinite line with the symmetric potential the 
amplitudes t and r are parametrized as follows {59! ]: 

t = cos6e iS , r = ism6e i5 } (2.59) 

where 9 and 5 are the real functions of u. Further we have 

det S(u) = t 2 (u) - r 2 (uj) = e 2lS{uj) (2.60a) 

= (2.60b) 

For the real u the function det S(co) is defined by formula (I2.60al) . while ( 12. 6 Ob j) gives its 
continuation to the whole complex plane u. 

In order to consider in this approach the initial asymmetric configuration (1|2|3), one 
has to proceed from the symmetric multilayered configuration (1|2|3|2|1), calculate the re- 
spective total reflection and transition amplitudes, and find detS(u). Transition to the 
initial asymmetric configuration should be done only at the final stage of the vacuum energy 
calculation (see Sec. IIII A[) by moving the interface 3|2 to infinity. 



III. CASIMIR ENERGY 



We define the Casimir energy as the renormalized zero point energy of the electromagnetic 
oscillations for configuration at hand 

£=~5> ? - (3-1) 

{9} 

Here the sum sign implies the summation over the discrete part of the spectrum and the 
integration over its continuous branch. 

The summing up of the discrete eigenfrequencies will be performed by making use of the 
frequency equations. The integration over the continuous part of the spectrum cannot be 
carried out by applying the 'naive' formula 

J ' codu, (3.2) 

because it does not 'feel' the physical nature of the problem under consideration, namely, 
this integral does not depend on the specific form of the differential operator and the corre- 
sponding boundary conditions [31]. Therefore, the summation over the continuous branch of 
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the spectrum should be accomplished, as usual, by making use of the spectral density shift. 
The rigorous mathematical theory of the scattering problem gives the following expression 
for this function 61 1 



Ap(u) = p(u) - p (u) = -—6(u) (3.3a) 

7r aw 

1 d , , „. . Id, t(u) 

\ndetS(u) = — : — In 4-4- (3-3b) 



2iri du x 2iri doo t(u>) 

where S(u) is the S"-matrix in the spectral problem at hand. In Eqs. ( 13. 3 j) p(u>) is the 
density of states for a given potential (or for given matching conditions in the case of the 
compound media) and poiui) is the spectral density in the respective free spectral problem 
(for the vanishing potential or for the homogeneous unbounded space). In the case of the 



one-dimensional scattering the calculation of the spectral density was considered in |59l. 162 1. 
It is worth noting here that already at the level of the spectral density (13. 3p the removing 
of infinity is carried out, namely, the contribution to the vacuum energy generated by the 
unbounded Minkowski space is subtracted here. 

Taking into account all this we can work out Eq. ( 13.11) in detail as follows: 



h . r°° <i 2 k . r°° 
E i 2a ) = o5Z / 7o~\2 ^2^n(cr,k,a) + uj Ap(a,u,k,a) du 



- (a oo) . (3.4) 



Here k is the wave vector along unbounded dimensions k = (k x ,k y ); u n (a,k,a) are the 
frequencies appertaining to the discrete spectrum; the function of the spectral density shift 
Ap(w) is given by Eqs. ( I3.3P ; uo is the lower bound of the continuous spectrum. Summation 
over a = TE, TM in (13.41) takes into account two polarizations of the electromagnetic field. 

The further treatment of Eq. (13.41) is aimed at deriving a unique integral representation 
for the contributions of both the discrete part and the continuous part of the spectrum, the 
integration being carried out over the imaginary frequencies u = i(. 

To be definite, we shall consider the most 'richest' electromagnetic spectrum possessing 
the surface modes, waveguide modes, and the continuous part (see, for example, the configu- 
rations I and IV in Sec. Ill B lj) . As was shown there, the frequencies of the surface modes lie 
in the interval < u sm < U-(k), the frequencies of the waveguide modes belong to the band 
uj-(k) < u sm < u) + (k), and for the frequencies u of the continuous part of the spectrum we 
have oj + {k) < uj (see Fig. [2]). 
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Im (i) 




FIG. 2: (color online) The contours C+ and C_ in the complex uj plane which are used when going 
on to the integration over the imaginary frequencies. For simplicity, the spectrum possessing one 
surface mode u sm , two waveguide modes uj wg , and the continuous part u > uj + is considered. The 
cuts starting at the points uj- and uj + are shown by heavy lines. 

A. Transition to imaginary frequencies 



In what follows, we have to go out in the complex frequency plane u. For this the one- 
valued branches of the integrand in Eq. ( 13.41) should be chosen. The square root dependence 



k a (u) = \ e a (oj)— - k 2 



a 



1,2 



(3.5) 



results in branch points. For the usually accepted functions e(u) there appear four such 
points on real axes: ±u)-(k) and ±u + (k) (0 < oj-{k) < cj + (k)) (see Sec. Ill B lj) . These are 
the same points that separate the different parts of the spectrum. 

We draw two cuts on the real axes u: the first cut connects the points — w_ and and 
the second one starts at the point u + and goes to infinity (see Fig. [2]). The first cut passes 
round the eigenfrequency of the surface mode u sm from below. It is not crucial and this cut 
can go round u sm from above as well. Further, we assume that ki(oj) and k 2 (oj) in Eq. (13. 5p 
acquire real positive values at the upper edges of the cuts. 

The eigenfrequencies uj n (a, k, a) are unknown in an explicit form, we have only the respec- 
tive frequency equations ( 12. 32ft and ( 12.33!) or m another form ( 12. 35ft - ( 12.37ft . This dictates, 
without choice, the making use of the argument principle theorem from the complex analysis 
when summing over n in (13.4]) . When choosing a contour for integration in the complex u 
plane one should bear in mind that the left-hand sides of the frequency equations D(a,u) 
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considered as the functions of the complex variable u have a 'good' behaviour in the upper 
half-plane u only 



D{a,u) = 1 - r y iak \ ck 2 = ^u 2 -u\. (3.6) 

This forbids from the very beginning the deformation of the contour in such a way that it 
includes the whole imaginary axes Imw. An admissible contour is obviously the following 
one: the integration should be performed along the circles C £ of the radius e with e — > 0, 
each circle surrounding a single root of the frequency equation (see Fig. [2]). 

When treating the contribution of the continuous part of the spectrum to the vacuum 
energy one should keep in mind the following. Upon substituting ( I3.3bl) in Ap(u)) the 
transition amplitude from (12. 54p 

«») = -o^r ' (3 - 7) 

in vacuum energy (13. 4 p there appear terms which are proportional to the distance between 
the slabs 2a. These contributions are generated by the factor e 2lk2a and they belong to 
the internal energy of the medium 2 in the configuration under consideration. These terms 
are removed from the Casimir energy by the subtraction indicated in Eq. ( 13.4[) . Hence, in 
calculation of the vacuum electromagnetic energy the denominator D(u) in Eq. (13. 7\\ 'works' 
alone. 

Now the vacuum energy (13 .4p acquires the form 



/ s hsr^ f°° kdk I 



d , \ r° , d , D(a,u) 

aoooo——[nD{a } ui)+ / duu— In— -. r 

c £ du J u+ du D(a,u) 



(3.8) 



Further, we represent the contribution of the continuous spectrum to the vacuum energy 
(the second term between the square brackets in Eq. (I3.8P ) in the form of the contour 
integrals as well. According to the Cauchy integral theorem, we can write 

duo uj-j— In D(a, tu) = 0, (3.9a) 
r du 

duu^-\nD(a,u) = 0, (3.9b) 

q_ dui 

where the contours C + and C_ are shown in Fig. [2j Here we have taken into account 
that the complex roots of the equation D(a,u) = lie only in the lower half-plane u and, 
respectively, the complex roots of the equation D(a, u) = are placed only in the upper 
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half-plane u. 6 Such a lay-out of the complex eigenfrequencies is in a direct relation with the 
choice of the time dependence in the solutions to the Maxwell equations in the form e~ wt 
(see Sec. ). In view of relations (13. 9p . we can express the contribution of the continuous part 
of the spectrum in ( J3.8|) as the integrals along the contours C + and C_ which have no section 
(cj + , oo). One can easily see from Fig. [2]that along the contours C £ the contributions of the 
discrete spectrum and continuous spectrum are mutually canceled. 7 Besides, at the length 
(0, u + ) the contributions of the continuous part of the spectrum generated by D(a, u>) and 
D(a, uj) are mutually canceled too. 

As a result, we arrive at a very compact and clear formula for the Casimir energy 

E(2a) = -J2 d(lnD(a,i(,k), o = TE, TM, (3.10) 

27r „ Jo 21T Jo 



where 



D(a, i(, k) = l- rl{i(, k) e 4aK2 ^ k \ a = TE, TM, 



«i(C, fc) + K2{(, k) £2«)«l(C) fc) + ei«)«2(C S *0 



«a(C> fc ) = **o(*C, fc ) = \/ e a [%Q-r + ^ 2 > « = 1, 2 . 



c 2 



In obtaining formula (I3.10p the integration by parts was carried out with respect to the vari- 
able (: u = i(. Obviously, this formula holds for an arbitrary spectrum of electromagnetic 
excitations sustained by the configuration under consideration. Let us remind that the sym- 
metric configuration is considered with nonmagnetic materials {ei{ui),fi\ = 1\e2{uj),H2 = 
l\ei(u),jj,i = 1). 

Differentiation of the vacuum energy (13.101) with respect to the gap width 2a gives im- 
mediately the Lifshitz formula for the Casimir force at zero temperature [l| 

F(2a) = ~l^af = -2^So kdk So dC K2iC ^ ? [r " (< ' k)eiaK2iC ' k) ~ " ' (3 ' 12) 
Here a = TE, TM as before and the same notation ( 13. lip is used. 
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A simple example of the complex frequencies in solutions to the Maxwell equations for an open system (a 
perfectly conducting sphere in vacuum) was discussed in Ref. 
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7 In earlier papers dealing with other configurations the analogous cancellation was brought out by making 



use of the complex fe-plane (see, for example, 63J, [6J]). However the Lifshitz formula is presented in 
terms of the variable iQ = u. 
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In order to reproduce exactly the Lifshitz formula (2.9) in Q| for a vacuum gap (62 = 
fi2 = 1) between the identical slabs a new variable p = p(k) should be introduced 



Taking into account that p(k = 0) = 1 and 



P 2 C 3 
dk = — — dp , 

c A 

we easily convert Eq. ( 13.1 2p to the form 

fe 2 poo poo 

F(2a) = - — ( 3 d( p 2 dpJ2K 2 «,P)e iap(/c -lV , (3.14) 
c Jo Ji CT 

where 

\V £ i«) - 1 +P -PJ \V £ i«) ~ l+P 2 -e 1 (i()p / 

IV. CONCLUSION 

Our derivation of the Lifshitz formula (13.101) is accomplished in the framework of the 
rigorous spectral summation method. We believe that in this approach the answers are 
given to many questions in this field which have been waiting for their solutions for a long 
time. First of all, it is simultaneous incorporation, on the same footing, of the discrete and 
continuous parts of electromagnetic spectrum, consistent treatment of the branch cuts in 
the oj plane, consideration of the complex eigenfrequency location and so on. The principal 
novelty of our approach is a special choice of appropriate passes in the contour integrals 
which are used for transition to imaginary frequencies. For the first time contours of this 



type were employed in our preceding paper 65[ dealing with a close problem. 8 The main 



idea of our approach is easy-to-grasp. It is sufficient to look attentively at Fig. [2] and the 
mechanism of mutual cancellations between contributions of the discrete and continuous 
parts of the spectrum becomes evident as well as the final result, Lifshitz formula (I3.10p . 
Without doubt, it is an important advantage of our consideration. 



By the way, it is these contours that had to be used in our preceding paper [66| considering dielectric 
cylinder without dispersion 
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By making use of the contours C + and C_ shown in Fig. [2]one can easily 'explain' in what 

3.101) and (13.111) may be 'derived' by summing up only the 



way the correct Lifshitz formulae ( 
discrete eigenfrequencies 

(flay ,llJ[ 



18j]). Let us ignore the branch cuts in the u plane, the 
complex eigenfrequencies (quasi- normal modes), and the bad behavior of the left-hand side 
of the frequency equation D(a,u) in the lower half-plane u. After that one may 'glue' the 
sections of the contours C + and C_ which are along the positive Re u axis and drop them. 
As a result, we obtain a unique contour involving the whole Imcu axis and a semi-circle of 
a large radius in the right half-plane. Using this contour in the first integral between the 
square brackets in Eq. (I3.8P we arrive immediately at the final Eq. (I3.10p . Proceeding in 
an analogous way one may 'derive' the Lifshitz formula in the scattering formalism alone 



mi2Qi). 

It should be noted here that the comparative analysis of the contributions to the vacuum 
energy generated by surface, waveguide, and photonic modes has recently been carried out 



by M. Bordag [67] . 



In the Lifshitz formula ( I3.10p written in terms of the integral over imaginary frequencies 
there are no traces of the electromagnetic spectrum details in the problem at hand; in other 
words, this formula has the same universal form for all physically admissible material media. 
Presumably, it implies that the problem of finding the vacuum energy with allowance for 
the material characteristics of media is, on its nature, a statistical problem. Therefore, 
an appropriate mathematical method here would be the dissipation-fluctuation theorem 
formulated, from the very beginning, in terms of imaginary time. One can anticipate that 
this may simplify the transition to imaginary frequencies in the Lifshitz approach. 

Many important topics are left beyond our consideration. One of them is the allowance 
for the dissipation when calculating the vacuum energy. There are many approaches to this 
problem (see, for example, |68f and more recent papers [69j]) but till now it is unclear how 
to come to an agreement with rigorous spectral summation method which is based on the 
real eigenfrequencies in the pertinent spectral problem. 
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